aarrdims=size(aarr)
slice=aarrdims[3]-1

print, 'contour plot set to the temporal slice', slice

slicemax=max(aarr[*,*,slice])
slicemin=min(aarr[*,*,slice])
nl=80
Lz=1.4
Lx=5.0
nz=240
nx=120
zarray=(Lz/nz)*indgen(nz)
xarray=(Lx/nx)*indgen(nx)

print, 'x dimension length', Lx
print, 'z dimension length', Lz

lev=grange(slicemin,slicemax,nl)
contour, aarr[*,*,slice], xarray, zarray,/iso, /fi, nl=nl
colorbar, range=[slicemin,slicemax], /horizontal, pos=[0.77,0.1,0.87,0.95], xtickformat='(F6.1)', xaxis=2, charsize=1.5

end
